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Abstract 

Assuming that neutrinos are Majorana particles, in a three generation framework, current and 
future neutrino oscillation experiments can determine six out of the nine parameters which fully 
describe the structure of the neutrino mass matrix. We try to clarify the interplay among the 
remaining parameters, the absolute neutrino mass scale and two CP violating Majorana phases, 
and how they can be accessed by future neutrinoless double beta (Qv/3/3) decay experiments, for the 
normal as well as for the inverted order of the neutrino mass spectrum. Assuming the oscillation 
parameters to be in the range presently allowed by atmospheric, solar, reactor and accelerator 
neutrino experiments, we quantitatively estimate the bounds on mo, the lightest neutrino mass, 
that can be inferred if the next generation 0vf3f3 decay experiments can probe the effective Majorana 
mass (m ee ) down to ~ 1 meV. In this context we conclude that in the case neutrinos are Majorana 
particles: (a) if rriQ > 300 meV, i.e., within the range directly attainable by future laboratory 
experiments as well as astrophysical observations, then m ee > 30 meV must be observed; (b) if 
mo < 300 meV, results from future 0^/3/3 decay experiments combined with stringent bounds on 
the neutrino oscillation parameters, specially the solar ones, will place much stronger limits on the 
allowed values of mo than these direct experiments. For instance, if a positive signal is observed 
around m ee = 10 meV, we estimate 3 <> m /meV <^ 65 at 95 % C.L.; on the other hand, if no 
signal is observed down to m ee = 10 meV, then mo ^ 55 meV at 95 % C.L. 

PACS numbers: 14.60.Pq, 13.15.+g, 14.60. St 



* Electronic address: nunokawa@ift.unesp.bi 
^Electronic address: teves@charme.if.usp.bi 
^Electronic address: zukanov@if.usp.br 



I. INTRODUCTION 



During the last few years, a significant amount of information on the size of the neutrino 
oscillation parameters have been gathered. Most of what we currently know about these 
parameters rely on evidences of neutrino flavor transformation that have been collected by 
the experimental observations of solar Q as well as of atmospheric neutrinos ||. The evi- 
dences coming from solar neutrino data have been strengthened by the recent neutral current 
measurement at Sudbury Neutrino Observatory (SNO) [§] whilst those coming from atmo- 
spheric neutrinos have also been strongly supported by the K2K accelerator based neutrino 
oscillation experiment |4[]. Furthermore, the negative results of reactor experiments || also 
impose stringent limits on some oscillation parameters. 

Assuming that only three active neutrinos participate in oscillations in nature, indepen- 
dent of whether neutrinos are Dirac or Majorana particles, current and future neutrino 
oscillation experiments can determine at the most six out of the nine parameters which 
completely describe the neutrino mass matrix, i.e., two mass squared differences (Amf 2 , 
Am^), three mixing angles (6 I 12 , 613, 623) and one CP violating phase (S) which parametrize 
the Maki-Nakagawa-Sakata (MNS) || leptonic mixing matrix. See, for instance, Refs. |7], |3| 
for recent discussions on the determination of these oscillation parameters. 

However, if neutrinos are of the Majorana type, there remain three non- oscillation param- 
eters, which can not be accessed by oscillation experiments. They are the absolute neutrino 
mass scale, which can be taken as the lightest neutrino mass, and two extra CP violating 
Majorana phases || |H], [DJ. It is well known that neutrinoless double beta (Qv(30) decay 
experiments can shed light on these non- oscillation parameters. 

The 0uf]f3 decay is a process which can occur if and only if neutrinos are Majorana 
particles [0. A positive signal of 0z//?/3 decay always implies a non- zero electron neutrino 
mass |13| even if it is not induced by the exchange of a light neutrino but by some other 
mechanism such as the one in supersymmetry models with broken i?-parity |14[]. In this 
work, we assume the simplest possibility to be true, that 0v(3f3 decay process is induced 
only by the exchange of a light neutrino. 

It has been abundantly discussed in the literature the relationship between the signals in 
Ovflfi decay experiments and oscillation phenomena, see for e.g., Ref. [ 15fl . So far, a large 
amount of effort has been made to constrain the oscillation parameters from the observation 
or non-observation of 0v(3f3 decay as well as to predict the possible range of the effective 
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Majorana mass in 0v(3(3 decay experiments, m ee , from the allowed range of oscillation pa- 
rameters fl5|l . 

In this paper, we take a different point of view. We examine how well we can constrain 
the three non-oscillation parameters by future 0vf3f3 decay experiments, considering that 
the oscillation parameters will be soon precisely determined (or constrained) by current and 
future oscillation experiments. We discuss the interplay among these parameters and the 
observable signal in future Qv(3f3 decay experiments for the normal as well as for the inverted 
ordering of the neutrino mass spectrum. In particular, presuming the oscillation parameters 
to be in the range presently allowed by the atmospheric, the solar and the reactor neutrino 
experiments, we examine what can be concluded about these parameters in the case of either 
a positive or a negative signal is obtained in future 0vf3f3 decay experiments fll6| . 

So far, no convincing signal of 0^/3/3 decay has been observed, rather only an upper bound 
on m ee 

m ee < 350 meV, (1) 

which comes from the result of the Heidelberg-Moscow Collaboration |T7J exists. Recently, 
an experimental indication of the occurrence of decay has been announced |18| but 

since this result seems to be controversial |19|], we do not discuss it in this work. 

There are many proposals for future Qvj3(3 decay experiments to go beyond the bound 
given in Eq. (@), those include GENIUS ||, CUORE 0, EXO MAJORANA || and 



NOON p4] . It is expected that in the initial phase of the proposed GENIUS experiment J2G 



the sensitivity to m ee can be as low as ~ 10 meV, going down to ~ 2 meV, if the 10 ton 
version of the experiment is implemented. In this work, we will try to be optimistic and 
consider that future experiments will eventually inspect m ee down to ~ 1 meV. 

The absolute neutrino mass scale is also independently constrained by the tritium decay 
experiments, which can directly measure the electron neutrino mass, obtaining the upper 
bound |25| 

m Ve < 2200 meV. (2) 
The proposed KATRIN experiment aims to stretch the current sensitivity down to ~ 340 



meV 26 . We also take this into consideration in our discussions. 



This paper is organized as follows. In Sec. |I|, we describe the theoretical framework on 



which we base our work. We first discuss in Sec. Ill, the dependence of the Qvj3{3 signal 



on the lightest neutrino mass, and second in Sec. IV we discuss how the dependence of m e 
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on mo is related to the two CP phases, a\ and a^. In Sec. [V] we discuss how m™ in , the 



minimum possible value of m ee , depends on m and #13. Finally, in Sec. [VI], we discuss 
how the upper as well as the lower bounds on m depend on the solar neutrino oscillation 



parameters. Sec. [VI I| is devoted to our discussions and conclusions. 



II. THE FORMALISM 



In this section we discuss the theoretical frame work we will rely upon in this work. 



A. Mixing and mass scheme 

We consider mixing among three neutrino flavors as, 







V\ 




= U 











(3) 



where v a {a = e, /1, r) and Vi{i = 1, 2, 3) are the weak and the mass eigenstates, respectively, 
and U is the MNS mixing matrix, which can be parametrized as, 

iS 



C12C13 



S12C13 



size 



-S12C23 - Ci2S23Sl3e i5 C12C23 - Si2S23Si3e i<5 S23C13 
Sl2'S23 — c 12 c 23 s 13 e *' 5 ~Cl2 s 23 — Sl2C23 s 13 e *' 5 C23C13 



(4) 



where s^- and Cy, correspond to sine and cosine of 9 it j. We define the neutrino mass-squared 
differences as Am^ = m| — m 2 , where Am^ = Am 2 2 is relevant for the solutions to the solar 
neutrino problem, and Am 2 tm = |Am 23 | — |Am 2 3 | is relevant for the atmospheric neutrino 
observations. 

Current atmospheric neutrino data indicate that 
Am 2 tm = |Am 2 3 |^(2-4)xl0- 3 eV 2 
which combined with nuclear reactor results [H] imply 

sin 2 13 < 0.02, 



sin 2 2^3 ^ 0.9 - 1, 



(5) 



(6) 



while the various solar neutrino experiment results strongly suggest that the so called large 
mixing angle (LMA) MSW solution with parameter in the range fl], |3|, |27| 

(2 - 40) x 10~ 5 eV 2 , tan 2 # 12 ~ 0.2 -0.8, (7) 



Am| 



l 12 



will prevail as the explanation to the solar neutrino problem. We will admit throughout this 
paper that the actual values of these parameters will be confirmed within the above ranges 
by future neutrino oscillation experiments. Besides further constraining the oscillation pa- 
rameters given in Eqs. it is expected that these experiments also can probe the CP 
violating phase S and determine the neutrino mass spectrum (sign of Am^) (7]. 

In this work, we denote the lightest neutrino mass by uiq. Then using Am| and Am^ tm 
as defined above, we can describe the two possible mass spectra as follows, 

(a) Normal mass ordering: 

mi = m Q , 



m 2 = Wmo + Am 2 



0) 



m 3 = V m o + Am| + Aml trn . (8) 



mi = Jml- Am 2 Q + Am ltm, 



(Jo) Inverted mass ordering: 



m 2 ee ^/m^ + Amf tm , 

m 3 == m . (9) 

In this manner Am\ 2 = Am 2 & for both mass ordering and Am^ = ±Am^ tm where the +(— ) 
sign indicates normal (inverted) mass ordering. In Fig. [I], a schematic picture of the mass 
ordering we consider here is shown. 



B. Effective Majorana mass and 

Assuming that the Qv(3f3 decay process occurs through the exchange of a light im v < 10 
MeV) neutrino, the theoretically expected half-life of the Ovfifl decay, Ty 2 , is given by |28 



[T^]- 1 = G°»\M 0v \ 2 ml e , (10) 



where G 0u denotes the exact calculable phase space integral, Mq u consist of the sum of the 



Gamow- Teller and the Fermi nuclear matrix elements defined as in Ref. |28j and m ee is the 
effective Majorana mass defined in Eq. (JTTf) below. 

It is known that the evaluation of the nuclear matrix elements suffers from a large uncer- 
tainty depending on the adopted method used in the calculations. As we can see in Table 



II of Ref. ||28|| , the evaluated half-lives for a given nucleus and a given value of m ee typically 
vary within a factor of ~ 10 comparing the largest and smallest predicted values. This im- 
plies a factor of ~ 3 difference between the minimum and the maximum values of m ee when 
extracting it from the results of Qvf3(3 decay experiments, which in fact directly measure or 
constrain not m ee but the value of Ty 2 . This is clear from Eq. (|10"D. In Sec. [V|, we will 
consider for our estimations a somewhat optimistic uncertainty of a factor ~ 2 instead of 3, 
assuming future improvements on the evaluation of the nuclear matrix elements. 
The effective Majorana mass, m ee , is given by, 

m ee = \m x XJl x + m 2 Ul 2 + m 3 U* 3 \, 

I 2 2 2ioi 2 2 2 2ia-\ I /1 1 \ 

= \mic 12 c{^e 1 +m 2 s 12 c l2i + m 3 s n e (11) 

where we have chosen to attach the CP violating phases to the first and third elements. 
Note that ot\ and a 3 must be understood as the relative phases of U e \ and U e3 with respect 
to that of U e2 . The range of these phases are 

< e*i < 7r, < a 3 < 7T. (12) 

As it is known, the value of m ee can be perceived as the norm of the sum of three vectors 
m®, m@) and m*fj in the complex plane whose absolute values are given by, 







= If/ 2 

— I°el 


mi 


2 2 
= 777,iC 12 c 135 








= If/ 2 

— 1^62 


m 2 


2 2 

— m 2 s 12 c 13 , 








= If/ 2 
— I^e3 


m 3 


= m 3 s 2 13 . 


(13) 



Explicitly, m ee is expressed as 

m L = [ m ie cos 2«i + + ""4e cos 2a 3 ] 2 + [m^ sin 2cti + sin 2a 3 ] 2 , 
= [m£>] a + [^i?] 2 + [mS? + 2 { m£>m£? cos2 ai 
+mj; 2) m^ ) cos2a 3 + m^mfj cos[2(«i - a 3 )] } , 

244, 244, 24,or 222 o 

= m 1 c 12 c 13 + rn 2 s l2 c l3 + m 3 s 13 + 2 { mim 2 c 12 s 12 c 13 cos2«i 

+m 2 m 3 s 12 c 13 s 13 cos2a 3 + mim 3 c 12 c 13 s 13 cos|2(ai — a 3 )J } . (14) 

We can clearly see from Eq. ([14] ) that m ee is invariant under the transformation 

(ai,a 3 ) -> (7r - «i,7r - a 3 ), (15) 
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which allows us to further restrict the range of (011,0:3), without loss of generality, to 

< «i < 7T, < o 3 < 7r/2. (16) 

We note that a± and/or 03 different from and tt/2 imply CP violation. 
The maximum possible value of m ee , denoted by m™ ax , is given by, 

= (miCi 2 + m 2 s 2 12 )c 2 13 + m 3 s 2 13 . (17) 

which occurs for ot\ = 03 = 0. On the other hand, the minimum possible value of m ee is 
zero only when the three vectors m|' (i = 1, 2, 3) can form a triangle. This can occur when 
the condition, 

™g<5>^ (* = 1,2,3), (18) 

is satisfied. When these three vectors can not form a closed triangle, which includes the 
case when one of them is null, the minimum value is given by twice the length of the largest 
vector minus the sum of the norm of all three vectors, 

m™ = 2max{m«,m( 2 e ),m( 3 e )} - [mW+mfl+mf']. (19) 

The values of (oti, 03) which lead to such a minimum are («i, 03) = (7r/2, 0), (tt/2, tt/2) or 
(0, tt/2) when respectively m^, mS 2 J or m^J is the largest contribution. 

C. Some useful extreme limits 

To help the comprehension of our discussions in the following sections, let us review here 
the approximated expressions for m™ ax and m™ m for the two extreme cases of the absolute 
neutrino mass scale: vanishing m and very large mo compared to Aml tm . The neutrino 
oscillation parameters are assumed to lie in the ranges given in Eqs. (fj)-(|7|)- 

(i) Vanishing mo limit: 

For the normal mass ordering we have, 

m- n , m?r - mf] =F m<$ * sj^s 2 ^ T (20) 
where the — (+) sign corresponds to m™ in (m™ ax ). 
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For the inverted mass ordering we have, 

m- in ~ mgj> - m% ~ J Arnicas 26 12 cj 3 , (21) 



w%T - m£> + m% ~ jAm^c^. (22) 



(ii) Large m (> y/Aml^) limit: 

For the normal as well as for the inverted mass ordering, 

<f - ™£ - mg> - rog ~ m (4 cos2fl 12 - s 2 13 ), (23) 
m™ ax ~ m$ + mf} + mf} ~ m . (24) 



III. DEPENDENCE ON THE LIGHTEST NEUTRINO MASS 

In this section we examine how the effective Majorana mass, m ee , depends on the lightest 
neutrino mass, mo, and clarify the importance of each of the individual contributions, m^, 
m@J and to m ee . We present in Fig. |2] for normal (upper panels) as well as inverted 
(lower panels) mass ordering, the maximum and minimum values of 77Zgg clS db function of m 
for vanishing #13 (left panels) and sin 2 #13 = 0.02 (right panels). In the same plots we also 
show the individual contributions of m^, mfj and by dashed, dotted and dash-dotted 
lines, respectively. 

Let us first discuss the case of normal mass ordering. As we can see from the plots in the 
upper panels of Fig. [|, for smaller values of m , rrv^J is the dominant contribution whereas 
for larger values of m , mfj dominates over the other contributions. For the typical values 
of the oscillation parameters allowed by the solar, atmospheric neutrino and reactor data, 
mf} is almost always the smallest contribution to m ee , being subdominant in m ee and it can 
only be important when there is a large cancelation between m^J and m@) . 

It is clear that a strong cancelation in m ee can occur if at least two of , m^J and m® 
are comparable in magnitude. Just by comparing their magnitudes in the plots in Fig. |2|, 
we can easily see for which values of mo a strong cancelation in m ee can occur. For the case 
if #13 = 0, m ee can be zero only at one particular value of m (see Fig. |2| (a)) whereas for 
the case if #13 7^ m ee can be zero for some range of m (see Fig. 0(b)). See Sec. for more 
detailed discussions about the dependence of m™ in on m and #13. 

In the case of inverted mass ordering, the situation changes significantly. Here m^J > 
m ee >> m ee an d m ™ n always must satisfy the condition, 

m™ n > v /Am2 tm cos20 12 ~lO meV, (25) 
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for any value of mo for current allowed parameters from solar and atmospheric neutrino data 
and no complete cancelation in m ee is expected as we can see clearly from the plots in the 
lower panels of Fig. [|. Therefore, if no positive signal of 0uf3/3 is observed down to ~ 10 
meV, inverted mass ordering can be excluded as long as neutrinos are Majorana particles. 



IV. DEPENDENCE ON THE LIGHTEST NEUTRINO MASS AND CP PHASES 



Let us next discuss how the m ee dependence on m is related to the two CP phases, 
ax and a^. We present in Fig. ^| iso-contour plots of m ee in the m -ai plane for vanishing 
#13, for the normal mass ordering, for some values of the solar parameters taken within the 
allowed ranges given in Eq. (|7]). Note that, in this particular case, there is no dependence 
on the phase a.3 in our choice of parametrization, which is clear from Eq. (|IT). 



Fist we note that in the plots, iso-contours are symmetric with respect to ai/it = 0.5. 
Second we note that if m ee is smaller than certain value, m^ rt , the iso-contours are closed 
which means that the possible values of both m and ot\ are bounded to some limited range, 
which does not include zero. The critical value of m ee under which the contour is a closed 
one is given by 

m^f ~ mfj = \jAm 2 Q s 2 12 ~ (1 - 10) meV, (26) 

this dependence is illustrated in Fig. |3|. If a positive 0u/3{3 signal is not observed down to 
these values, this will imply either that mo as well as the CP phase a\ are bounded to 
the limited range within the closed contours shown in Fig. [3] or that neutrinos are of not 
Majorana type particles. 

We show in Fig. |]the same information as in Fig. |3](b) but for sin 2 #13 = 0.01 (left panels) 
and 0.02 (right panels) for three different values of a% = 0, tt/4 and 7r/2 in the upper, middle 
and lower panels, respectively. As we can see from these plots the qualitative behaviors of 
the contours are very similar to those in Fig. ||](b). This is due to the fact that m^J , which is 
the only term that carries the contributions of sin 2 #13 and 03, is subdominant compared to 
the other two elements in m ee . The effect of a non-zero a% is to cause some displacement of 
the position of the symmetry line of the plots from around ati/ir — 0.5 to somewhat smaller 
values. 

Let us here mention the case where mo can be independently measured by another exper- 



iment such as the KATRIN tritium decay one. In this case, it is possible to constrain 



the CP phase ot\ by comparing the measured values of m ee and mo provided that m > 340 
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meV, the maximum sensitivity of KATRIN. If a Ou/3/3 decay experiment measures m ee sig- 
nificantly smaller than mo measured by KATRIN, this would imply a non-zero CP phase 
ax. This is because for the m values relevant for KATRIN, m ee ~ m if ol\ ~ 0, but m ee 
can be as small as ~ 0.1 x m if o>\ ~ 7r/2 and if the largest allowed value of #12 from the 
current LMA allowed region is realized (see Eqs. ([23|) and (|24D ). However, it would still be 
difficult to say something definite about the value of ot\ for 10 < m /meV < 340. 



Our plots in Figs. |3] and f| are in agreement with the conclusion presented in Ref. p9 
that is, either a positive or a negative results in a Ovflfl decay experiment can constrain 
m and «i but, since the possible values of «i will always include tt/2, a non-zero value of 
a\ cannot be interpreted as an evidence of CP violation, even if a positive signal of 0^/3/3 
decay is observed. Unfortunately, to be able to say anything more definite on the CP phase, 
independent precise information on mo is unavoidable. Moreover, nothing can be concluded 
about the value of a^. 

We show in Fig. |5| the same information as in Fig. ||] for the inverted mass ordering. Since 
there is no significant dependence on # 13 or on a% in this case, we show only the curves 
for vanishing # 13 . We note from these plots that there is no lower bound for m as long 
as m ee < 50 meV, and moreover, ot\ is less constrained if compared to the case of normal 
ordering, independently of the values of the other neutrino oscillation parameters. 

For the case of the normal mass ordering, we have also investigated if the uncertainties in 
the determination of the solar parameters as well as on m ee can wash out the determination 
of a non-zero CP phase ati, expected by the closed contours in Figs. [3] and [|. For four different 



central values of m ee , assuming 30% uncertainty in their determination (see Sec. [VI| for a 
detailed explanation), we have obtained the region in the (Am 2 2 , tan 2 6*12) plane where a\ 
can be constrained to a non-zero value for (i) sin 2 #13 = and (ii) sin 2 #13 = 0.02. This is 
shown in Fig. [| where we also have assumed for each point in the (Amf 2 , tan 2 #12) plane 
a 10 % uncertainty in the determination of these two parameters. In this plot we have 
indicated by crosses the set of solar parameters used in Fig. |](a)-(d). 

We observe that for the case of #13 = (see Fig. |6| (i)), determination of a non-zero 
CP phase ot\ is possible as long as one can reach the sensitivity m ee < 5 meV for the 
paramerter set (d) and m ee < 1 meV for sets (b) and (c). On the other hand, in the case 
where sin 2 # 13 = 0.02, as we can see from Fig. || (ii), a better sensitivity in m ee is required 
to establish non-zero a\ value for the same parameter set. This is because when #13 is non- 
zero, third element m^J, which contains another CP phase, a^, whose value is assumed to 
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be unknown, comes into play and this could wash out more efficiently the determination of 
non-zero a± when compared to the case where #13 = 0. 



V. DEPENDENCE ON THE LIGHTEST NEUTRINO MASS AND 9 13 

In this and the next sections we focus on the relation among m ee and some of the yet 
undetermined mixing parameters. Let us start by discussing the dependence of m™ m on m 
and #13. Here we discuss only the case of normal ordering since the dependence on # 13 for 
the inverted case is quite small. We present in Fig. [!], the iso-contour plots of m™ m in the 
s 2 3 -mo plane for some different choices of the oscillation parameters. 

For vanishing #13, one can have m™ m = if m£) = m^J, i.e. , mic 2 2 = rri2sl 2 , which is 
equivalent to, 



m = - n 12 yAml ~ 3 meV, (27) 
' cos 26*12 1 



where we have computed the numerical estimate using the best fitted values of the param- 
eters from the latest solar neutrino data (see Fig. |2|(a) and Fig. 0(b)). Within the current 
allowed range given in Eq. (|7|), the possible values of mo for which we can expect strong 
cancelation are in the range m — (0.9 — 12) meV, as we can confirm with the plots in Fig. [7]. 

For non-zero values of #13, as we can see clearly from Fig. |7|, m™ in can take zero for some 
range of m . We can also see that there is a critical value of #13 for which m™ m is zero with 
vanishing m . Such a value of # 13 can be easily estimated by solving m 2 c 2 3 ,s 2 2 = m 3 sf 3 with 
itlq — > which is equivalent to, 



A 13 



Ami 



2 

© „2 



^atm 



0.03, (28) 



\ Am 2 

for the best fitted parameters from solar as well as atmospheric neutrino data. We note that 
this is close to the current upper bound on s 2 3 allowed by the CHOOZ result 0. Letting 
Arrig, Am 2 tm and s\ 2 take any value in the region consistent with the solar and atmospheric 
neutrino observations given in Eq. (|5]) and (0), the range of s\ 3 for which strong cancelation 
can occur is s\ 3 ~ 0.01 — 0.20, which is again consistent with our results in Fig. [F[ 

We observe that if s 2 3 is smaller than the critical value given in Eq. (|28| ) , the value of mo 
can be strongly constrained to some limited range around the value of mo given in Eq. (ETj) . 



provided that future Oz//3/3 experiment can probe a m ee value as small as ~ yAmgS^ 
independent of whether a positive or negative signal of Ou/3/3 is observed. 
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VI. CONSTRAINING m USING SOLAR NEUTRINO DATA 



Finally, let us discuss how we can constrain mo using the solar neutrino parameters. Let 
us first analyze the case where a positive signal of 0^/3/3 is observed. The value of m ee 
has to be extracted from the experimental measured half-life, T®Y 2 , of the decaying parent 
nucleus by comparison with the theoretical prediction which rely on nuclear matrix element 
calculations. This means that the experimental value of m ee has to be expressed as an 
interval obtained using the maximum and minimum matrix element predictions. 



As mentioned in Sec. II B, typically, there is a factor 3 difference among the results of 



the matrix element evaluations according to different model assumptions p8[ . In order 
to take such large theoretical uncertainty into account in our estimations, we will first 
assume that the experimental measured value of m ee will be extracted using the mean 
between the minimum and maximum values of the matrix element calculations, then attach 
30 % uncertainty around this value, which correspond to a factor ~ 2 between the smallest 
and the largest theoretically allowed m ee values for a given value of the observed half-life. 
This is a somewhat optimistic but reasonable assumption. Hopefully improvements on the 
understanding of the underlining nuclear physics effects can decrease this uncertainty even 
further. 

In Fig. H we show the iso-contours of upper (m^) and lower (m™ 111 ) bounds on mo 
in units of meV in the tan 2 #i 2 -Am 2 2 plane for the case where a positive signal of 0^/3/3 
is observed with central values m ee = 10, 5, 3 and 1 meV with a 30 % uncertainty. In 
these plots we have used Am| 3 = 3 x 10~ 3 eV 2 and sin 2 #13 = 0. We do not present plots for 
m ee > 10 meV, since in these cases the upper and lower bound can be analytically estimated 
as we will see below. 



When mo > yAm 2 tm the upper bound on mo does not depend much on Am|, but 
essentially only on # 12 (see Eq. (^) in Sec. [11 C| ), it is given as, 

~ (29) 

and the lower bound m™ 111 ~ m ee , independent of the solar parameters in the region com- 
patible with the LMA MSW solution to the solar neutrino problem. 

We observe that, as can be seen in Fig. that as m ee decreases, the upper bound lines, 
which mainly depend on tan 2 #12, are shifted to larger values of 612, decreasing m™ ax for a 
given set of (tan 2 #i 2 , Am 2 2 ). The lower bound lines, on the other hand, depend more on 
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Amf 2 and there are some regions where no lower bound is obtained. For m ee > 10 meV or 
m ee < 1 meV there is always a lower bound found inside the currently allowed LMA MSW 
region, whereas for 1 < m ee /meV < 10 this is not true. The appearance of these no lower 
bound bands comes from the fact the in these regions the solar mass scale alone is sufficient 
to explain the positive signal observed, even for vanishing tuq. 

In Fig. |]we repeat the same exercise but for sin 2 #13 = 0.02. The most significant effect 
of a non zero sin 2 9 is is the increase of the size of the no lower bound band, which can in 
some cases stretch over the entire LMA allowed region, otherwise the behavior of the upper 
and lower bound lines is qualitatively similar to the previous case. 

Let us next consider the case where no positive signal is observed. In Fig. [K] we present 
the same information as in Fig. |8| but for the case where no positive signal of Ou/3/3 is observed 
down to m ee = 10 meV (a), 5 meV (b), 3 meV (c) and 1 meV (d). It is assumed that when 
no positive signal is observed, for a given bound on the half-life time, the bounds on m ee 
are extracted using the smallest nuclear matrix element prediction which lead to the largest 
m ee value (see Eq. (0)). We can see that qualitative behavior of the iso-contours for the 
upper bound is similar to that in Fig. |8| but it is different for the trend of the lower bound 
curves. Compared to the case where positive Ovfifi signal is obtained, it is more difficult to 
put lower bound on tuq. This can be easily understood from Fig. § (a). We note that unless 
we can constrain m ee down to ~ 3 meV, no lower bound on ttlq is obtained. In Fig. 11 we 



present the same information as in Fig. pi] but for sin 2 # 13 = 0.02. Again, the qualitative 
behavior of the iso-contours is similar to the case of sin 2 9 13 = 0. The difference from the 



previous case presented in Fig. |10| is that the constraint on itlq become somewhat weaker, 
leading to larger upper bounds and smaller lower bounds for a given set of (tan 2 #12, Amf 2 ). 

Finally, let us also comment the case of inverted ordering. For this case, we can easily 
estimate the upper as well as lower bound from the analytic expressions as well as from the 
lower panels of Fig. |2|. Let us look at Fig. ^ (c) and (d). First of all, if the inverted ordering 
is the case, the observed value of m ee in 0uf3(3 decay experiment must be larger than the 



value given in Eq. (|25|), as already mentioned in Sec. pLTT| . If the observed m ee is smaller than 



Am 2 tm c 2 3 (see Eq. (]22D ) then there is no lower bound on m whereas the upper bound is 
given by the same expression for the case of normal ordering, by Eq. (p9[). If the observed 



m ee value is larger than yAm 2 tm c 2 3 , then the upper as well as lower bounds are given by 
the same expression as in the case of normal ordering, which is already discussed in this 
section. 



13 



VII. DISCUSSIONS AND CONCLUSIONS 



We have studied how the non- oscillation neutrino parameters which can not be extracted 
from the oscillation analysis, the absolute neutrino mass scale and two CP violating Ma- 
jorana phases, can be accessed by the positive or negative signal of future Qv/3P decay ex- 
periments. We have carried out this by choosing various different set of mixing parameters 
which are varied within the parameter region currently allowed by the solar, atmospheric as 
well as by reactor neutrino experiments. 

In the future, the KATRIN experiment expects to directly inspect m down to ~ 340 
meV ||26|| , while there are several proposed astrophysical measurements on the temperature 
perturbation in the early universe imprinted in the cosmic microwave background radiation, 
such as the ones that can be performed by MAP (Microwave Anisotropy Probe) |31| and 
Planck 1 32], which can probe mo down to ~ 300 meV |33], where the expected sensitivity 
suffers from the uncertainty coming from cosmological parameters. It is expected that future 
supernova neutrino measurements can probe m down to at most ~ (2-3) eV ||34|| . 

Our analysis permit us to conclude that if these experiments measure m > 300 meV, 
then either a positive signal of 0vf3f3 decay compatible with m ee > 30 meV must be observed 
in the near future or neutrinos are Dirac particles. This is simply because m™ m ~ mo cos 26±2 
(see Eq. (|23|) in Sec. [II Q) and cos2#i2 ^0.1 from the current allowed LMA parameter region 
given in Eq. (|7|). 

On the other hand, if these experiments do not observe any positive signal of mo ^ 300 
meV, results from future Ovflfl decay experiments combined with more precise values of 
the neutrino oscillation parameters, specially the solar ones (Am 2 2 , tan 2 #12), which can be 
precisely determined by the KamLAND experiment ||, will place more stringent bounds on 
mo for the Majorana case provided that the sensitivity of m ee < 30 meV is achieved. To 
be more specific, if a positive signal is observed around m ee = 10 meV (assuming a 30% 
uncertainty on the determination of m ee ), we estimate 3 < m /meV < 65 at 95 % C.L.; 
on the other hand, if no signal is observed down to m ee = 10 meV, then m < 55 meV 
at 95 % OL. Allowing for a more optimistic sensitivity, a positive signal observed around 
m ee = 3 meV or no signal seen down to 3 meV would mean mo < 25 meV at 95 % C.L. 
These bounds can be improved by a better determination of tan 2 6*12, Am 2 2 and sin 2 6*13 as 
can be clearly seen in Figs. |8f- |TT], as well as by the reduction of the uncertainties in the 
theoretical calculations of the nuclear matrix elements. 
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We finally conclude that it is possible to constrain the CP violating phase ai\ to values 
around ir/2 if: (1) mo is large enough to be detected by KATRIN or by astrophysical 
observations and future Oz//3/3 decay experiments observe m ee close to its minimum value or 
(2) future 0u(3j3 decay experiments can achieve a sensitivity on m ee < 5 meV, depending on 
the values of the solar parameters and on the uncertainty on m ee (Fig. ||), independently of 
whether a positive or a negative signal is observed. Unfortunately, nothing can be known 
about «3. 

Despite of the fact that the parameter region inspected in this work differs from the one 



considered in Ref. [59|, we have found, in agreement with this reference, that evidence for 
CP violation cannot be observed by future OV/3/3 decay experiments, since the possible values 
of cxi, whenever it can be constrained to some non-zero value, always include 7r/2, unless mo 
can be independently determined by some other experiment. 
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FIG. 1: Mass ordering considered in this work. 
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FIG. 2: Maximum and minimum possible values of 777-gg ctS cl function of mo indicated by the 
thick solid curves for sin 2 #13 = (left panels) and 0.02 (right panels) for normal (upper panels) as 
well as for inverted (lower panels) mass ordering. We have fixed the other mixing parameters as 



Am 2 2 



of 771 



5 x 10~ 5 eV 2 , tan 2 6»i2 



0.35 and [Am^l = 3 x 10 3 eV 2 . The individual contributions 



ee ; m ee 



and rrv e J are also shown by dashed, dotted and dash-dotted curves, respectively. 
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FIG. 3: Iso-contour plots of m ee for sin 2 #13 = for the normal mass ordering in the a\ — rriQ 
plane. We have fixed the other relevant mixing parameters as (At7i 2 2 , tan 2 #12) = (2xl0 -5 eV 2 , 
0.2) in (a), (5xl0~ 5 eV 2 , 0.35) in (b), (5xl0~ 5 eV 2 , 0.5) in (c) and (4xl0~ 4 eV 2 , 0.6) in (d). 
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sm 2 13 = 0.01 sin 2 13 = 0.02 
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FIG. 4: Same as Fig. ||(b) but for sin 2 9\z = 0.01 (left panels) and 0.02 (right panels), for 
a3 = (upper panels), 7r/4 (middle panels) and ir/2 (lower panels). The dashed vertical lines 
mark ax/w = 0.5. 
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FIG. 5: Same as Fig. ||| but for the inverted mass ordering. Here we have fixed the atmospheric 
mass scale to Am^ = — 5 x 10~ 3 eV 2 in (a), to Ara^ = — 3 x eV 2 in (b) and (c) and to 
Am| 3 = -1.3 x 1(T 3 eV 2 in (d). 
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FIG. 6: Region in the (Amf 2 , tan 2 # 12 ) plane where ot\ can be constrained to a non-zero value, 
indicated by an arrow, for given central values of m ee = 1,2,3 meV, taking into account 30% 
uncertainty in the determination of m ee and 10 % uncertainty in A?n 2 2 as well as in tan 2 6»i2 • The 
set of solar parameters used in Fig. H(a)-(d) are indicated by crosses. The allowed region for the 
LMA MSW solution at 1, 2 and 3 a are shown by the shaded area (adopted from Ref. p0| ). 
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FIG. 7: Iso-contour plots of m™ m in the s^ 3 — mo plane for different choices of the mixing 
parameters for the normal mass ordering. We note that in the region delimited by thin solid curves 



mf e m = 0. 
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FIG. 8: Iso-contour plots of upper (m™, solid curves) and lower (m™ 11 , dashed curves) bounds 
of mo in units of meV in the tan 2 612 — Am 2 2 plane for the case where a positive signal of 0v(3(3 
is observed with central values, m ee = 10 meV (a), 5 meV (b), 3 meV (c) and 1 meV (d). We 
assume 30% uncertainty in the determination of m ee . We fix the other mixing parameters as 
Am| = 3 x 10 -3 eV 2 and sin 2 #13 = 0. The allowed region for the LMA MSW solution is the same 
as in Fig. |6|. 
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FIG. 9: Same as Fig. § but for sin 2 9 i3 = 0.02. 
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FIG. 10: Same as Fig. |8| but for the case where no positive signal of Qv(3f3 is observed down to 
m ee = 10 meV (a), 5 meV (b), 3 meV (c) and 1 meV (d). 
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FIG. 11: Same as Fig. 10 but for sin^ 13 = 0.02. 
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